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Abstract 

We discuss the analysis of polarization experiments with particular 
emphasis on those that measure the Stokes parameters on a ring on 
the sky. We discuss the ability of these experiments to separate the 
E and B contributions to the polarization signal. The experiment 
being developed at Wisconsin university is studied in detail, it will be 
sensitive to both Stokes parameters and will concentrate on large scale 
polarization, scanning a 43° degree ring. We will also consider another 
example, an experiment that measures one of the Stokes parameters 
in a 1° ring. We find that the small ring experiment will be able to 
detect cosmological polarization for some models consistent with the 
current temperature anisotropy data, for reasonable integration times. 
In most cosmological models large scale polarization is too small to be 
detected by the Wisconsin experiment, but because both Q and U are 
measured, separate constraints can be set on E and B polarization. 

PACS numbers: 98.80.-k, 98.70.Vc, 98.80.Es 



Email address: matiasz@arcturus.mit.edu 



1 Introduction 



Temperature anisotropies in the cosmic microwave background (CMB) can 
be considered one of the best probes of the early universe. The higher pre- 
cision measurements that will become available in the near future could po- 
tentially lead to a precise determination of a large number of cosmological 
parameters |1|, 0] . But the temperature anisotropies are not the only source 
of information, Thomson scattering of the CMB photons will lead to a small 
degree of linear polarization at a level detectable by the future generation of 
experiments. 

The detection of polarization will provide information about our universe 
that is not available in the temperature data. For example it can provide an 
accurate determination of the ionization history of the universe [^. Gravi- 
tational waves and vector modes also leave a specific signature in the polar- 
ization ^. The correlation function of the Stokes parameters Q and U 
can be used to test the causal structure of our Universe and thus provides a 
direct test of inflation [^. The polarization power spectra produced by the 
current family of topological defect models are also very different from that 
produced by inflationary models 0. 

The two future satellite missions, MAP and PLANCK, will measure po- 
larization as well as temperature anisotropies. The extra information in 
polarization will help improve the constraints these missions can put on cos- 
mological parameters For most parameters the expected error bars de- 
crease by a factor of two and for those related to the reionization history of 
the universe the accuracy can be several times better. 

Several ground based polarization experiments are now also under way. 
Some of these experiments plan to measure the Stokes parameters in pixels 
on a ring around the north celestial pole, which simplifies the pointing of the 
telescope. The properties of the experiment being developed at Wisconsin 
University are discussed in detail. This experiment will concentrate on 
large angular scale polarization measuring both Q and f/ in a 43° ring. We 
also consider another example, an experiment that measures only one Stokes 
parameter in a 1° ring. 

Recent theoretical developments 0, ^ have shown that rather than de- 
scribing polarization in terms of Q and U it is more natural to introduce two 
scalar quantities E and B. The correlation functions of this variables is what 
theories most easily predict and the treatment in terms of E and B takes full 
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account of the spin 2 nature of polarization. 

The ring hke polarization experiments appear as an ideal testing ground 
for this new formalism because they are being performed right now and also 
because their simple geometry will enable us to study them in full detail 
and build intuition on this new polarization formalism. We will study for 
example how this kind of experiments can separate the two components of 
polarization E and B. Our treatment will also permit us to explore how 
changes of the experimental set up, like a change of the ring size, affect the 
experiments ability to detect cosmological polarization. 

The outline of the paper is the following: in §2 we review previous results 
on polarization and describe the way in which general polarization experi- 
ments can be analyzed. In §3 we discuss specific properties of the ring like 
experiments and in §4 we study parameter estimation issues in the context 
of this kind of experiments. We make a summary and discuss our results in 
55. 



2 Correlation Functions 

The CMB radiation field is described by a 2 x 2 intensity tensor lij . The 
Stokes parameters Q and U are defined as Q = (In — /22)/4 and U = /12/2, 
while the temperature anisotropy is given by T = (In + /22)/4. The fourth 
Stokes parameter V that describes circular polarization is not necessary in 
standard cosmological models because it cannot be generated through the 
process of Thomson scattering. While the temperature is a scalar quantity 
Q and U are not. They depend on the direction of observation n and on the 
two axis (ei, 62) perpendicular to n used to define them. If for a given n the 
axes (ei, 62) are rotated by an angle ip such that e'l = cosip ii + sinip 62 and 
62 = — sin ip ei + cos ip 62 the Stokes parameters change as 

Q' = cos2iIj Q + sm2ilj U 

U' = -sm2^ Q + cos2ip U (1) 

To analize the CMB temperature on the sky it is natural to expand it 
in spherical harmonics. These are not appropriate for polarization, because 
the two combinations Q ±iU are quantities of spin ±2 [T^. They should be 
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expanded in spin- weighted harmonics ±2^™ 0, 

Tin) = J2^T,im^imin) 

Im 

{Q + iU){n) = ^ a2,im 2yim{n) 

Im 

{Q-iU){h) = Y.<'-Vm-2YUn). (2) 

Im 

To perform this expansion, Q and U in equation (§) are measured rela- 
tive to (61,62) = (661,6^), the unit vectors of the spherical coordinate sys- 
tem. There is an equivalent expansion using tensors on the sphere The 
coefficients ±2aim are observable on the sky and their power spectra can 
be predicted for different cosmological models. Instead of ±2aim it is con- 
venient to use their linear combinations aE,im = —{cL2,im + o_2,/m)/2 and 
ciB,im = — (02,/m — a-2,/m)/2i, which transform differently under parity. Four 
power spectra are needed to characterize fluctuations in a gaussian theory, 
the autocorrelation between T, E and B and the cross correlation of E and 
T. Because of parity considerations the cross-correlations between B and the 
other quantities vanish and one is left with 

('^T,im'^£^,«m) — ^m,m'Cch (3) 

where X stands for T, E or B, {■ ■ ■) means ensemble average and Sij is the 
Kronecker delta. 

For the purpose of this work it is more useful to rewrite equation (|^) as 

Im 

Q{n) = — dE^lmXi^lm + iaB,lmX2,lm 
Im 

U{n) = — ^ 0'B,lmXi^lm — iaE,lm-^2,lm (4) 
Im 

where we have introduced Xi^im{n) = ( 2'i^m + -2^m)/2 and X2,im(n) = 
( 2Yim - -2Yim)/'2. They satisfy Xli^ = Xi^i^.^ and Xli^ = ~X2,i-m which 
together with aE,im = and a^./m = (i*B,i-m make Q and U real. 
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In fact Xi^imin) and X2^imi'h) have the form, Xi^im{n) = y(2/ + l)/47r 

Fi,imiO) e'^^ and X2,/™(n) = ^{21 + l)/47r F2,z^(^) e^"^^ F(i,2),«™(^) can be 
calculated in terms of Legendre polynomials |Q 



(/-2)!(/-m)! /-m^ 1 _ , 

+ {l + m)^^Pr^A^ose)] 
sm (7 



{I - 2)\{l - m)\ m 
(/ + 2)!(/ + m)!sin2 0' 
+ (/ + m)P,"^i(cos^^)]. 



-(/ - l)coseP;'"(cos( 



(5) 



Note that -F2,«m(^) = if m = 0, as it must to make the Stokes parameters 
real. 

The correlation functions can be calculated using equations (01) and (HI), 
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where 1 and 2 stand for the two directions in the sky ni and n2. These 
expressions can be further simplified using the addition theorem for the spin 
harmonics [O, 

E s^Yi*^{ni) s2Yim{n2) -- 



' 21 + 1 
An 



S2 



(7) 



where P is the angle between hi and n2, and {4'i,4'2) are the angles {ee^e^) 
at hi and n2 need to be rotated to become aligned with the great circle going 

subroutine that calculates this functions is available at 



http : / / arct urus . mit . cdu / mat iasz / CMBFAST 
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through both points. In the case of the temperature equation(|^) gives the 
usual relation, 

2/ 4- 1 

{T,T,) = J2^CTiPi{cosf3). (8) 
I 47r 

For polarization the addition relations for Xi^im and X2jm are calculated 
from equation (J^), 

2/ + 1 

J2 X*, ,„(l)Xi,z™(2) = — [Fi^dP) cos 2^1 cos 2^2 - F2,i2(/3) sin 2^i sin 2V'2] 

m 

2/ + 1 

5]X2%(l)X2,/™(2) = ^—[Fi,;2(/3) sin 2^1 sin 2^-2 -F2,z2(/?) COS 2^1 COS 2V^2] 

m 

J2XliJl)X2,lrni2) = Z^^[Fi,;2(/5) siu 2^^! COS 2^2 + F^^P) COS 2^^! siu 2^2] 
m 

2/ + 1 

5: ol^:;(l)Xi,z^(2) = ^^Fi,,o(/3) cos 2.^2 

9/ + 1 

^ o>^:;.(l)Xi,,„(2) = -^^_Fi,,o(/5)sin2^2 (9) 



where we can equivalently write -Fi,io(/5) = y - 2)!/(/ + 2)!P;^(/5). 

The correlations in equation (1) with the sums given by equation (|^) 
are all what is needed to analize any given experiment. These relations are 
simple to understand, as pointed out in the natural coordinate system 
to express the correlations is one in which both (61,62) reference frames are 
chosen to be aligned with the great circle connecting the two directions (1 
and 2); in that case we have |^ 

I 47r 

{Ur{l)Ur{2)) = J2'^^[CbiF,,12{P) - CeiF^MP)] 

I 47r 

(T(l)a(2)) = E^C'c.^i,o(/3) 
I 47r 

(T(l)f/,(2)) = (10) 

the subscript r here indicate that the Stokes parameters are measured in 
this particular coordinate system. We can use the transformation laws in 
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equation to write {Q,U) in terms of {Qr,Ur) and then using equation 
( p!0D for their correlations one can recover our final result (given by equations 
® and (I)). 

When analyzing an experiment we can arrange the measured values of T, 
Q and U in a. vector x = (Ti, Qi, Ui, ...Ti, Qi, Ui...), the subscript labels the 
pixel. The measured Stokes parameters at each pixel have two contributions, 
the first coming from the cosmological signal we are interested in measuring, 
given by equation (^) and the second from the noise in the detectors. The 
correlation matrix of x is (xjXj) = Cij = S^- + Njj, where the signal cor- 
relation matrix S is given in equations (|^) and (^) and N is the correlation 
matrix of the noise. With this correlation matrix the full likelihood can be 
calculated, 

L{^\Cxi) oc exp[-ix^C-ix], (11) 

here Cxi stands for the complete set of power spectra that describe the theory 
under consideration, Cxi = {Cti, Cei, Cbi, Cci} ■ If the set of power spectra 
are given in terms of a set of parameters, the likelihood can be maximized 
to find parameters that best fit the data. 

3 Ring Experiments 

There are several ground based polarization experiments now under way. 
Some of these experiments plan to measure the Stokes parameters Q and U 
in pixels that form a ring on the sky. The simple geometry of the sky patch 
is perfect for understanding the relation between E and B polarization and 
the Stokes parameters. It is the aim of this section to analyze this kind of 
experiments. 

Rather than doing the analysis in terms of correlations in real space as was 
suggested in the previous section, in this case it is simple to diagonalize the 
signal correlation matrix. If the noise is uncorrelated and equal from pixel to 
pixel (ie. the noise correlation matrix is proportional to the identity matrix) 
we simultaneously diagonalize the noise and signal correlation matrices. This 
will allow us to make a more detailed analysis of the characteristics of this 
type of experiments. 

Let us assume that the experiment measures the Stokes parameters in 
a number (Npi^) of pixels on a ring of radius 9 with a gaussian beam of 
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width Ofuihm- We will take the noise correlation matrix to be diagonal, 
{^{Q,u)iN(Q,u)j) = o'pSij and uncorrelated between Q and U, N(^Q^ijy are the 
noise contribution to the measured Q and U in pixel j. If = 2n{j — 1)/Npix 
denotes the angle along the ring, the Stokes parameters in each pixel are given 
by (equation (|) plus the noise contribution) 



Im 



Qj 

Uj = -J2^{21 + l)/A7rBi^[aB,imFi,imi0) 

Im 

-taEMF2MidW'^^' + Nuj, (12) 

here Bi^^ encodes the information the beam and scan pattern of the exper- 
iment. In an appendix we make the derivation of these functions for the 
Wisconsin experiment. In what follows we will take Bf^ = exp [— /(/ + l)o"fe] 
with (Tft = ^/«,/im/2\/21n2 which is a good first approximation. 

Diagonalizing the signal correlation matrix is very simple and intuitive: 
if one considers the Fourier transform of the data then each mode will pick 
a particular value of m in equation (|T^). In particular let us consider the 
transformed data set 



Y: ^{21 + l)/A7,Bi,[aE,iA{e) + iaB,ikFl{e)] + N, 

l>\k\ 

1 ^ ^ " -ik<f>j 



k 



TV • ^ 

^ "ptx j 



= V(2/ + l)/4nBi,[aB,ikFl,{e) - taE,ikFUe)] + (13) 

l>\k\ 

where k runs from —Npix/2 < k < Npix/2. To get to the last expression we 
have used that, 

i(m—k)<f)j ^ ^ ^ ^i{m—k)2'iT /Npix^j 



_ g*{''"-fc)27r^^^]^ _ ^i{m-k)2iT/Npi^^ 
-^pix^m,k+nNpix 

(14) 
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where n is any integer. The sum in equation (|14|) is only different from zero 
if = m or if it differs by a multiple of Npi^. To obtain equation (^) we 
ignored the terms that "leak" power from higher harmonics, only k = m 
is considered. This is a reasonable assumption because higher harmonics 
are suppressed by beam smearing. In fact we should always choose Npi^ 
large enough so as to oversample the beam and not loose the information 
in the smaller scales, this immediately guaranties that the aliased power is 
negligible. 

In equation (|13[) we have labeled the Fourier transforms of the noise con- 



tributions N^Q jj-^ which under the assumption of uniform uncorrelated noise 
satisfies, 

(iVfiV^') = 

(iV^*iVj) = (15) 

where wj} = ap/Npi^;. 

We will take the data set to be the Fourier coefficients rather than the 
Stokes parameters. The signal part of each Fourier component receives con- 
tributions only from multipoles with m = k, equation then implies that 
only components with the same value of k will be correlated. The correlation 
matrix is block diagonal, 

l>k 

{u'*u''') = 6,,,,Y.i^i + ^)/^^BUCE,iFl,ie) + Cp,iFl,,ie^ 

l>k 
+Sk,k'Wp^ 

{U'^^Q^') = 5k,k'Y.^{2l + l)/A-nBl{CE,i + CB,i)Fi,iu{e)F2,ik{e) 

l>k 

{Q''*U^') = -{U^*Q^'), (16) 

"*" means complex conjugate. An important point is that this transforma- 
tion does not rely on the particular shape or other property of the power 
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spectra, if the CMB is a statistically isotropic and homogeneous random 
field taking the Fourier transform will always make the correlation matrix 
block diagonal. The correlation matrix is hermitian, and the cross correla- 
tion between IJ^ and is imaginary. It will be more convenient to change 
variables and use iV^ IJ^ ^ the data set will then be = {Q^,ilJ^). The 
correlation matrix for each k is given by equation (|1^) but without the i in 
the cross term. 

In the Fourier domain the likelihood becomes a product of independent 
gaussians, one for each fc, 

L{ii\Cxi) oc n exp[-^x^CfcXfc], (17) 

YvdetCfc 2 

where the matrices are 2x2 given by equation (|I6|) and f means trans- 
pose and complex conjugate. Clearly if only one of the Stokes parameters is 
measured then the likelihood is just the product of several one dimensional 
gaussians and if the temperature is also included in the analysis turns 
into a 3 X 3 matrix. 

Equation ([T6|) shows that if a ring experiment only measures one of the 
Stokes parameters, either or tj'^ it will not be able to separate the E 
and B contributions in a model independent way because both contributions 
enter in the expressions summed together. This is similar to what happens if 
one is interested in separating the contributions from density perturbations 
and gravitational waves to the temperature anisotropics. It is only when 
one assumes a shape for the power spectra that one can fit the observed 
temperature spectra as a sum of these two separate contributions and infer 
the presence of gravity waves. In a similar way if we only measure and 
we assume different shapes for both Cei and Cbi we could determine each 
contribution. On the other hand if both Stokes parameters are measured then 
we could separate the E and B contribution directly as they enter differently 
in each correlation (multiplied by Flij^{6) and F2n.{6) in Q and the other 
way around in U). 

Another interesting property of equation (|TBp is that the = Q mode 
does not receive any contribution from the B channel because -F2,/o(^) = 0, 
the converse is true for U . For example a non zero signal in f/° implies 
the presence of B and thus of gravity waves or vector modes. Of course 
this determination suffers from a large cosmic variance, because we are only 
measuring one realization. 
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3.1 Small scale experiment 

Let us now analize how the signal in each mode depends on the power spec- 
trum. We will do this first for a small ring experiment that measures only Q 
on a ring with 9 = 1°, 6fu,hm = 0.22°. We selected a small ring scanned with 
a fraction of a degree resolution to make the experiment sensitive to small 
scale polarization. We will choose Wp^ = (0.9/ii^)^, corresponding to a re- 
ceiver noise of 1.2mK/\/Hz and three weeks of observations. This numbers 
are chosen to be representative of what a ground experiment might achieve. 
The correlation matrix is 

(g'^^g^') = 4 (2^ + l)/^7T{Wi,lkCE,l + W2,lkCB,l) + Sk,k'Wp\ (18) 

Both signal and noise parts are diagonal and we have introduced the window 
functions Wi^ik = Bf^ F^^ki^) and W2,ik = Bff^ F^^iki^^) which describe how 
each Fourier mode depends on the underlying power spectra. Let us consider 
what the experiment would measure if the underlying cosmological model was 
standard CDM. We will assume there are no gravity waves or vector modes 
so we will only consider E type polarization. We will generalize our analysis 
later. 

Figure |1| shows the window functions Wi^ik for several values of k. Only 
multipoles I > \k\ contribute to Fourier component k. As |fc| grows the peaks 
of the window functions move to higher / (we only show results for the A; > 
modes but they also apply to the A; < ones). The window functions are 
very broad in / space, because their width is inversely proportional to the 
smallest dimension of the sky patch one is observing. For large values of / 
the window functions are cut off by the effect of the beam smearing, Bfj.. 

Figure |l| also shows that the most important contribution to the signal in 
this experiment comes from the second polarization acoustic peak at / ~ 350, 
the first peak at / ~ 100 contributes to a smaller extent. Higher peaks are 
inaccessible because of beam smearing. 

An interesting point to note is that the k = 2 window function peaks at 
low / (the lowest / in the polarization expansion is by definition / = 2), so it is 
this mode rather than /c = or 1 that gets the contribution from the largest 
angular scale perturbations. To understand why let us consider a small patch 
of the sky containing the 6 = 1° ring and a polarization field that is constant 
over the patch. This polarization field is produced by the very large angular 
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Figure 1: Window functions for the small Ving experiment for several values 
of k. The SCDM spectrum is also shown for comparison, we actually plot 
[l + \)CeiI'^t^ which is the relevant measure of power for a linear / scale (the 
normalization is arbitrary). 



scale perturbations, the low I modes. Without loss of generality let us call tx 
the direction of the polarization vectors, will be perpendicular to that. The 
Stokes parameters are measured relative to the (e^, basis which rotates as 
one moves along the ring. Although the polarization amplitude P is constant 
over the whole patch and the direction is always the Stokes parameters 
along the ring (measured in the (e^, e,/,) basis) are given by 

= P cos 2(f) j 
Uj = Psm2(pj, (19) 

and thus the k = ±2 Fourier modes will get all the signal. The coordinate 
system used to describe the Stokes parameters, which is rotating as one moves 
along the ring, and their spin 2 nature, which tells us how Q and U transform 
under this rotation, makes the large angular scale perturbations contribute 
only to the k = ±2 mode. This is another illustration of the spin 2 nature 
of polarization. 

Another interesting thing to consider is the effect of the size of the ring 
on the window functions. The window functions are oscillating functions of 
/ and changes in the ring size modify the period of these oscillations. In 
particular the separation between peaks scales as 1/6, and at the same time 
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Figure 2: Window functions for the small ring experiment for k = 4,8 and 
two ring sizes 6 = 1°, 0.8°. 

as the peaks become more separated they also become broader. This is most 
easily understood in the small scale limit. If we can neglect the curvature of 
the sky, which is an excellent approximation for this experiment, we can use 
Fourier transforms instead of spherical harmonics and write ^ 

Q(9) = (27r)-2 J dHE{l) cos 2(0 - 0,) e''^™^(*-<^') (20) 

where 6 labels the Fourier mode with [l^ + ily) = /exp (z0;). We have again 
considered only E type polarization and here {E(\i)E{l2)) = CeiiS^{12 — li), 
5^(x) is the Dirac 5— function. A straightforward calculation leads to 

Q' = -^J dHE{l)e-''^^ \[Jk+2{u) + Jk-2{u)], (21) 

where u = 19 and Jn{u) are Bessel functions. The correlation matrix is then 
given by 

(Q'*Q'') = h,k' J ^ Cei \[Jk+2{u) + Jk-2{U)]' (22) 

which means that the window function in the small scale limit is Wi^ik = 
l[Jk+2{u) + Jk-2{u)]'^. The effect of changing the size of the ring is straight 
forward to understand in this expression, each Wijk is only a function of m = 
W so changes in 6 just cause a stretching of the window function. A similar 
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Figure 3: Signal to noise in the different modes for the small ring experiment 
assuming w]^^ = {0.9 fiK)^ and COBE normalized SCDM as the underlying 
model. Two different ring sizes are plotted, in each case the effect of changing 
the beam width from 0.22°(upper curve) to 0.2° (lower curve) is also shown. 



analysis for the other set of window function yields W2,ik = \[Jki+2{u) — 
Jki-2{u)]'^- Beam smearing only adds a Bff. to the expressions of both W(^i^2),ik 
. Figure ^ shows two examples of these window functions for two ring sizes, 
the effect of scaling together with the Bfj^ factor is clear. 

In figure § we show the S/N for each of the Fourier modes, that is the 
ratio of the signal to the noise part of the correlation in (equation ITBI), as a 
function of k. As expected the signal to noise decreases with increasing k 
because of beam smearing. The effect of decreasing the ring size is evident 
in the figure, because the window functions for successive k become more 
separated and at the same time broader the signal in the experiment gets 
concentrated in fewer Fourier modes. Also note that there is a peak at /c ~ 7 
when the ring is ^ = 1° because this window function just matches the second 
acoustic peak and is then able to pick more power. The position of this peak 
is shifted to lower k when the size of the ring is decreased. 

It is important to recognize that most modes have a. S/N smaller than 
one, which will make the polarization hard to detect. Other models different 
from SCDM would produce a higher signal, but reducing the noise in the 
experiment would crucially improve its sensitivity to cosmological polariza- 
tion. On the other hand the overall S/N of the experiment not just mode by 
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Figure 4: Window functions for the Wisconsin experiment for several values 
of k. For comparison in the middle pannel we show the E spectra, (/ + 
l)C'_Ei/2vr, for SCDM but with = 1 (the normalization is arbitrary). 



mode is better, opening the possibility that cosmological polarization will be 
detected in the near future with this type of experiments. 

3.2 Wisconsin experiment 

Let us now consider the Wisconsin experiment, because it will actually mea- 
sure both Q and U on the ring it is in principle capable of distinguishing 
between E and B polarizations. The ring will be 6' = 43°, scanning the sky 
on the vertical at Wisconsin. It will have Ofwhm = 7° and Wp^ = {OAfiKY, 
the experiment is thus sensitive to large scale polarization. Unfortunately 
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the large scale polarization is very small in most cosmological models for the 
sensitivity of this experiment. 

Figure ^ shows some of the window functions for the experiment with 
Wi,,, = Bl Fl,{e), W,^iu = Bl and Ws,ik = Bf, F.^^e) F,,M. 

Note that the window functions for the cross correlation between Q and U 
can actually be negative. The properties of these window functions are the 
same as those discaused above for the small scale experiment, but note that 
because both the ring and the beam are much bigger now, the experiment is 
only sensitive to / < 25. 

We will illustrate the features of this experiment by considering different 
underlying toy models, for simplicity we will take the power per logarithmic 
interval of E and B to be constant, /(/ + 1)C(^e,b),i/'^'^ = ^'{e,b)- Figure ^ 
shows the Q and U S/N, the ratio of the signal to the noise contributions 
to the autocorrelations, for a model with a"^ = 2fiK'^ and no B. For com- 
parison we also show what is expected for SCDM with an optical depth due 
to reionization Tri = 1.0. Note that even for this large optical depth the 
expected signal to noise is very small. The smaller k modes have smaller 
S/N in the = 1 model than in the flat spectra model because large scale 
polarization power spectra decreases rapidly with / in all realistic models. In 
figure § we also show the E spectra for a model with = 1. All models 
where the universe reionizes at an early enough epoch have a peak in the 
large angular scale polarization spectra 0. For this model Ipeak ~ 20 and 
in general Ip^ak ~ 2^5^ with Zri the redshift at which the universe reionizes. 
We can see that the window function of the experiment are more sensitive to 
larger angular scale polarization and thus although the peak of the = 1 
model has an amplitude of ~ 6fiK^ the fiat spectrum model with amplitude 
2fiK^ produces a higher S/N for most k. 

The experiment can only detect cosmological polarization for extreme 
models with a very early reionization. On the other hand it has the virtue of 
being able to put separate constraints on E and B because it measures both 
Q and U. This is the feature of the experiment that we are most interested in. 
We will analyze simplified models with constant power spectrum for both E 
and B. This will allow us to study what kind of constraints a ring experiment 
can put on E and B separately, not relying on the different shapes of these 
two contributions. The approximation also makes visualization of the results 
easier and it is also straight forward to get results for any other noise level. 
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Figure 5: S/N for the Stokes parameters in the Wisconsin experiment. At 
low k for both Q and U the lower curve corresponds to COBE normalized 
SCDM with Tri = 1.0, the upper curves are for a model with cr|; = 2fiK'^ and 
al = OfiK^ 

one just needs to rescale the values of cf'^e,b) accordingly. 

4 Estimating Parameters 

A useful way of characterizing the amount of information a given experiment 
can provide is through the Fisher matrix. If one is interested in constraining 
a set of parameters s, the Fisher matrix is defined to be 

A, ^ C-f . (23) 

For example the minimum error bars that can be obtained on a given pa- 



rameter is Asj = yFjj"^. 

The Fisher matrix can be used to construct minimum variance quadratic 



estimators for the different parameters [p!3|, P 
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These estimators also have the virtue of being unbiased, {sj) = Sj. The 
estimator depends on the parameters one intends to estimate through the 
correlation matrix so some iterative approach may be needed to implement 
this method. This procedure can be seen as an iterative way of maximizing 
the likelihood fl^. On the other hand a reasonable choice of underlying spec- 
trum usually leads to a good estimator which has a slightly higher variance 
than could be achieved if the "correct" underlying model was used to build 
the estimator. Perhaps a more important problem is that one also makes a 
mistake in the error bars of the estimator which could lead to errors in the 
subsequent determination of cosmological parameters from the data. We will 
discuss this further in the following section. 

4.1 Small scale experiment 

Let us start with the simplest application of these formulas. Let us imagine 
that we are only interested in determining one parameter from the small ring 
experiment, we will assume that the shape of the spectrum is given and that 
we want to obtain an amplitude. We will parametrize the power spectrum as 
Cei = PCei, Cei is a fiducial spectrum which might be for example SCDM. 
The correlation matrix and its derivative are then, 

C{k,k') = {(3ak + w-^)6k,k' 
dC{k, k') 



dp 



akSk,k'- (25) 



where we have defined = J2i>k{'^^ + ^)/4:iTWi^ikCEi- Note that ak is both 
a function of the experimental set up, ie. the ring and beam size and of the 
fiducial cosmological model. 

As we are only estimating one parameter the Fisher matrix is just a 
number (F), 

F=IT. ^o f -i,r (26) 
2 V {(Jak + Wp y 

We can use this Fisher "matrix" to obtain the expected error bars on (3, 

A(3_ _ 1 
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Figure 6: Relative error in (3 for three different underlying cosmological mod- 
els as a function of the size of the ring 6. The two curves in each case are for 
beam size 0.22° on the top and 0.2° on the bottom. 



^Ekl/il + wp'//3a,y ^^^^ 

This can be interpreted as saying that the effective number of modes con- 
tributing informations is N^jf = 1^^1/(1 + Wp^ /Pak), where clearly the 
modes with a higher signal to noise contribute more information {S/N\k = 



We can use equation (|27|) to investigate the effect of changing the size 
of the ring 6, which enters which enters in the calculation of a^. Figure 
P shows A/3//3 as a function of 6 for three cosmological models. The E 
polarization spectra of these models is shown in figure ^ we choose standard 
CDM (SCDM) and two other models that fit well the available temperature 
data ||15|. One is an open model with = 0.85, fi^/i^ = 0.026, h = 
0.4 and n = 0.91. The other is a flat CDM model with fi^, = 0.05 and 
h = 0.3. It is interesting to note that while this two models have similar 
temperature spectrum up to the first acoustic peak, within the accuracy of 
the current CMB measurements, the amplitude of the polarization spectra 
differ significantly. Figure ^ shows that this experiment may be able to detect 
polarization in the foreseeable future, at least for some cosmological models. 
Increasing the sensitivity by reducing the noise and the angular resolution, 
thus going after the larger small scale polarization, will increase the prospects 
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Figure 7: E polarization power spectra for the three models discussed in the 
text. 



of the experiment as well. 

The quadratic estimator for beta also takes a very simple form, 
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T,k<^l/{Pak + Wp 



1^2 



(28) 



which is nothing more than an inverse variance weighing of each mode. It 
is straight forward to check that the ensemble average of /3 is actually P 
((|xfcp) = Pak + Wp^), it is an unbiased estimator. 

In this very simple example we can understand what happens if we con- 
struct the estimator using the "wrong" power spectra. Let us assume that 
the underlying model is SCDM model (ie. that the underlying (3 = 1 so that 
(|xfcp) = ak + Wp^) but that we use a different /5 to construct the estimator. 
The first thing that is easy to check is that {j3) = 1 always, the estima- 
tor is always unbiased even if we use the incorrect /3 to construct it. What 
happens is that the variance of the estimator increases if use the incorrect 
fiducial model, the variance in the estimator is given by 



J2k c^lKPak + Wpfx {ak + Wp )V(/5«fc + Wp 



1\2 



[Ek^l/ipak + Wp 



1^212 



(29) 



This function is shown in figure H, the variance has a minimum at (3 = 1. If 
one chooses an incorrect value of (3 the estimator is still unbiased but is not 
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Figure 8: Variance of the f3 estimator if the underlying model has {3 = 1 
but a fiducial model with a different (3 is used to construct the quadratic 
estimator. The curves have been scaled using the minimum variance value, 
when the /? = 1 model is used to construct the estimator. 

optimal, one could get more information out of the data (ie. reduce the error 
bars on j3). Note that the if /? ~ 2 which means we are overestimating the 
power in the underlying model by a factor of 2, the variance only increases 
by a few percent. Even for /3 ~ 4 the increase is less than 10 %. The curves 
in figure ^ depend on the signal to noise ratio of the experiment, for example 
the open model curve grows more slowly than the other two. This model has 
less power so more of the variance in the estimator comes from the detector 
noise and so it is less sensitive to the fiducial model. 

Unfortunately life is not so simple, it is not true that all the models we are 
interested in testing have a spectrum of the form Cei = 13Cei as is obvious 
from figure 0. The open and flat models were chosen to fit the existing tem- 
perature data [0 , this means that their first acoustic peak is approximately 
in the same place and is approximately of the same height. The acoustic 
oscillations in the photon-baryon fluid are responsible for both temperature 
and polarization peaks so if the temperature peaks of two models coincide in 
position, the polarization ones will also coincide. They not necessarily have 
the same amplitude, as can be seen in figure 0. To get the best results in this 
type of analysis we should choose Cei = I^Cei with Cei corresponding to one 
of the models that fit the existing temperature data. Perhaps polarization 
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Figure 9: Comparison between the underlying spectrum and the "recovered" 
one. In the top pannel we show C|;'"* and {I3)C°^^^. The bottom pannel 
shows the S/N for each k mode for the flat model and the scaled open one. 

can even help distinguish between some of these models before smaller scale 
temperature measurements do. 

In fact if we construct the quadratic estimator under the assumption 
that the underlying model is the open one, the mean value of 13 if the real 
model is the fiat one is {(3) = 3.033, figure || compares /(/ + l)C^^i* 1271 and 
{(5)l{l + l)C°^i^ 1271. The agreement is very good especially for the S/N curves 
for the two models, which is what the experiment is ultimately sensitive to. 

In any case (3 should be taken just as an indication of the power in po- 
larization, and because this experiment is sensitive to such a wide range in 
I it is better to assume some realistic spectrum shape given that the low 
/ polarization is so suppressed. If we do not do so and just assume a flat 
spectrum the weighing of the k modes in j3 will be very bad, leading to an 
estimator with a very large variance. Of course one could divide the / range 
in many bins so that the fiat approximation in each bin is reasonable, but 
the S/N in the experiment does not permit that yet. The best way to go is 
to assume a shape for the spectrum corresponding to some model that fits 
the existing temperature data well. 
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4.2 Wisconsin Experiment 



Let us now consider the Wisconsin experiment, with the assumption that 
both E and B spectra are constant, the correlation matrix for each k is 



where a^j = I]/>fc(2/ + 1) / iirWj^ik for j = 1,2 or 3 and 1 is the identity 
matrix. 

Note that there is a symmetry between {E,B) and {Q,U), that is to 
say interchanging E and B and Q and U leaves equation (|30| ) unchanged. 
Consider for example the previous signal to noise plot for a model with 
a% = 2^K'^ and no B (figure In the "conjugate" model having a% = 2fiK'^ 
and no E, the Q and U signal to noise plots would just be interchanged. The 
signal to noise plots for both parameters are different so at least in principle 
the E and B contributions can be separated. 

We have two parameters so the Fisher matrix is 2 x 2. we used it to 
compute the expected error bars on a"^ and in a grid of underlying models. 
Contour plots of the error bars are shown in figure 0. One interesting 
feature of these plots is that because of the symmetry discussed above both 
figures are essentially the same, one can go from one to the other by simply 
interchanging the axis. 

The other feature that deserves attention is that for a fixed value of one of 
the amplitudes, for example a"^ as one increases the other one (a^) the error 
bar on a"^ grows. What happens is that the determination of both amplitudes 
is highly correlated, and always the combination that is better constrained 
tends to be more aligned with the variable with the highest amplitude. We 
show the error ellipses in the (t|; — a% plane for two models in figure 0. The 
first model has a"^ = 2fiK^ and = 2fiK^, the second one a'^ = 3fiK^ 
and cr|; = IfiK^. Clearly because of the E — B symmetry in the correlation 
matrix if the underlying model has a"^ = the sum of both amplitudes (ie. 



(XfcX[,) 





(30) 
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Figure 10: Contour plots of the cr^^ error bars that could be obtained by 
the Wisconsin experiment for a grid of underlying models with E and B 
amplitudes between O^K^ and lO^K^. 



the total power) is much better constrained than their difference. When we 
consider a model with more power in the E channel the ellipses rotate so 



that (Tp is better constrained than a? 



In figure 12 we show the error bars expected for under the assumption 
that (T^ = 0, that is models that lay on the x-axis of figure |10[ We show 
two curves, one correspond to trying to determine both a'^ and a% and the 
other assumes we know there is no B polarization. This is actually a good 
approximation for most models, only defect models that usually have a large 
vector contribution have a large B contribution. These models do not seem 
to be in very good accord with the existing CMB and LSS data [jl6| . 



This experiment would be able to detect the cosmological polarization 
in an underlying model with a"^ ~ IfiK'^. A model with r^j = 1 could be 
detected, but this model has most of its large angular scale power outside the 
band where the experiment is most sensitive. The "reionization peak" in this 
case occurs at Ipeak ~ 20 with an amplitude GfiK'^, thus a smaller beam would 
greatly increase the sensitivity of the experiment to the reionization history 
of the universe, making the experiment better "tunned" for the position of 
the peak in / space. 

It is also clear that a ring experiment can put independent constraints on 
both E and B separately if it measures Q and U. Although the assumption 
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Figure 11: Contour plots of 68% and 95% confidence in the a% — a% plane 
that could be obtained with the Wisconsin experiment. The axis have been 
scaled using the amplitudes of the underlying models. Two different models 
are considered. 

that the spectrum is flat across all the band is not very good in this case, if 
the noise in the experiment was reduced then several bins in / could be used 
and our analysis would become more realistic. 

5 Summary and Discussion 

We have studied in detail polarization experiments that measure the Stokes 
parameters on a ring in the sky. We have shown that the Fourier transform 
of the Stokes parameters on the ring provides natural variables to analyze 
the results of the experiment, and have studied how the correlation matrix 
of this quantities depend on the E and B power spectra. We calculated the 
window functions and discussed how it depends on the size of the ring. 

Experiments that measure both Stokes parameters are capable of deter- 
mining separately the contributions from both E and B. They can do so 
without assuming a specific shape for each of the spectra because E and B 
contribute differently to Q and U autocorrelations. In fact for A; = the Q 
mode only depends on E while the U mode is only sensitive to B. 

We have calculated the Fisher matrices for these experiments and used 
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Figure 12: Expected error bars on a"^ for models with no 5 as a function of 
the underlying models a]^. The lower curve results if we assume there is no 
B rather than fitting for its amplitude. 



them to estimate their ability to detect the amplitude of cosmological polar- 
ization. We also constructed the minimum variance quadratic estimator for 
the parameters of interest and studied their properties in detail 

An experiment targeted at small scale polarization measuring the Stokes 
parameters in a small ring with a fraction of a degree resolution, may de- 
tect cosmological polarization in the near future. In fact, with sufficient 
integration it may help distinguish between models that have very similar 
temperature power spectra in the / range measured so far and so could add 
useful information. 

The Wisconsin experiment is mostly sensitive to large scale polarization 
but because the signal in this part of the spectrum is very small it will be 
harder for it to detect polarization. On the other hand the measurement of 
both Q and U will allow the experiment to put separate constraints on E 
and B. 

The prospect of detecting polarization is increasing quickly and in the 
near future we may have a positive detection. We will then be able to study 
very interesting characteristics of our universe, like its ionization history or 
the presence of vector and tensor modes. 
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6 Appendix 

In this appendix we will consider the window functions for the Wisconsin 
Experiment. This experiment will be pointing in the vertical direction scan- 
ning a ring in the sky with 6 = 43°. The instrument will be rotated in steps 
of 45° to measure both Q and U 0. The signal at each pixel is given by 

1 rto+St/2 f 

Xj = — dt dn'B{n{t),h')X{n) (31) 

Ot Jto-5t/2 J 

where X stands for the Stokes parameter being measured, j represents the 
pixel and St is the time interval the instrument is measuring each particular 
Stokes parameter before being rotate by 45° degrees. The beam is represented 
by the convolution with B{{n{t),n') and n{t) is the instantaneous position 
of the beam, in this case given by 6' = 43° and = ujst with ujs = 27r/day. 

We can use the expansion in harmonics (equation ^ to calculate the 
convolution with the beam. For a gaussian beam this gives the usual factor 
exp [— /(/ + l)cr5/2]. For example 

Qj = ^/rft[-E«i^,'-e-'('+^)'^'/2Xi,„(0,0(t)) 

Im 



+taBMe-'^'^'^''''^'X,,iUO, 4>m (32) 

and a similar expression applies to U. 

The time integration is simple, the dependence is only in the exp [im0(t)] 
factors inside X(i^2),«m- The window function then acquires another factor 
from 

1 f j^jm6(t) im<A, sin(mA0/2) 



dtJ'^^'' mA0/2 ^^^^ 



^This is one of the proposed strategies, a continuous rotation of the instrument is also 
being considered. 
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with A0 = ujs6t, the size of the pixel. The combined window is Bim = 
exp [— /(/ + l)(T^/2]sin(mA0/2)/(mA0/2). The importance of this last term 
depends on the size of the pixel, and can be made irrelevant by choosing 
small enough pixels. We ignored this factor in the main text for simplicity, 
but it can trivially be accounted for. 
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